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Summary 

The propagation of spinning mode sound was con- 
sidered for a cylindrical duct with sheared steady flow. 
The calculations concentrated on the determination of 
the wall optimum acoustic impedance and the maximum 
possible attenuation. Both the least attenuated and 
higher radial modes for spinning lobe patterns were 
considered. A parametric study was conducted over a 
wide range of Mach numbers, spinning lobe numbers, 
sound frequency, and boundary layer thickness, A cor- 
relation equation was developed from theoretical con- 
siderations starting with the thin boundary layer approx- 
imation of Everaman. This correlation agrees well with 
the more exact calculations for inlets and provides a 
single boundary layer refraction parameter which deter- 
mines the change in optimum wall impedance due to re- 
fraction effects. 

Introduction 

Noise attenuation data using inlet suppressors, with 
wall treatment only, have pointed out two phenomena 
which have required additional study. First, the noise 
attenuation was much larger than had been anticipated 
from flow duct tests or from axisymmetric sound propa- 
gation theory. Second, large attenuations were ob- 
served with much higher porosity (low acoustic resis- 
tance) acoustic liners than would have been anticipated 
from previous optimization studies. A recent paramet- 
ric study of the optimum wall impedance for the least 
attenuated spinning mode 1 '^ has shown that spinning 
modes could be anticipated to attenuate much faster than 
axisymmetric modes in an inlet with wall treatment 
only. This has helped explain the first point raised 
above. It was suspected that the bulk of the answer to 
the second point above involved the effect of the steady 
flow inlet boundary layer on the optimum wall imped- 
ance. This study was thus initiated to extend the work 
of Ref. 1 (spinning mode propagation with acoustic 
treatment) to include the boundary layer effects. This 
paper wall concentrate on the optimum w r all impedance 
and maximum possible sound attenuation and the corre- 
lation of these results with the pertinent- variables in- 
volved. 

The effect of boundary layers on sound propagation 
in acoustically lined ducts has been considered in the 
past for axisymmetric mode (nonspinning) propagation 
(e.g. , Refs. 2 and 3), and an extensive bibliography 


and discussion of the recent papers on this subject is 
given in Ref. 4. Spinning mode sound propagation with 
an axially uniform boundary layer was presented in 
Ref. 5. Comparisons of theory and experimental data 
were presented showing reasonable agreement. The 
above studies are very important since they demonstrate 
the necessity of considering boundary layer effects, but 
they do not include a parametric study of the problem 
which can be immediately used for suppressor design 
purposes. Current interest in sound propagation with 
sheared flow is apparent in two recent papers. ^ ^ 
Reference 6 presents a parametric study of the optimum 
impedance using the least attenuated axisymmetric mode. 
Some calculations are presented for spinning modes. 
Reference 7 presents spinning mode calculations but no 
optimization is attempted. 

The purpose of this paper is to present a parametric 
study and a correlation of the results for the effect of 
axially uniform boundary layers on the propagation of 
spinning mode sound in an acoustically lined cylindrical 
duct. A wide variation of steady flow Mach number, 
boundary layer thickness, sound frequency parameter, 
and spinning lobe number were considered. The em- 
phasis was on the least attenuated radial mode for the 
various spinning lobe numbers; however, the second and 
fifth radial modes for a seven lobe pattern were also 
analyzed. 

The analysis consisted of locating the optimum wall 
impedance (at which maximum possible sound attenuation 
occurs) for the previously discussed wide range of vari- 
ables. 

A major objective of this study was to obtain a cor- 
relating equation to unify the results of the parametric 
study. Since the correlation is quite complicated, ap- 
proximate expressions are developed that more clearly 
demonstrate the effect of the variables involved and that 
are sufficiently accurate for first estimates in acoustic 
liner design. 

Symbols 

A specific acoustic admittance (1/C) 

Aq specific acoustic admittance - uniform or slip flow 

B 0 optimum resistance coefficient for circular duct 
with slip flow (see Eq. (32) and Table I) 



optimum reactance coefficient for circular duct 
with slip flow (see Eq. (33) and Table I) 

c speed of sound, m / sec 

D circular duct diameter, m 

F boundary layer refraction function (see Eq, (27)) 
f frequency, Hz 

G boundary layer shape function (G = (y/5) ' ' for 
l/7^ power profile) 

1^ integral function across boundary layer (see Eq. 
<A-6)) 

I 2 integral function across boundary layer (see Eq. 
(A- 7) or Eq. (28) for l/7^ power profile) 

i V-i 

J Bessel function of first kind of order m 

m 

K factor in axial wave number, see Eqs. (3), (9), 
and (10) 

k a?/ c, m“ ^ 

M axial steady flow Mach number (function of radial 
coordinate) 

M q axial steady flow Mach number - free stream 
uniform value ‘ 

m spinning mode lobe number 

2 

P acoustic pressure, N/ni 

p part of acoustic pressure which is function of 
radial coordinate (see Eq. (3)), N/m 2 

Q 1 + iM^cr + ir) or 1 - MKq 
B amplitude of eigenvalue a 

<R ratio of optimum resistance to optimum reactance 
for uniform flow (no boundary layer), - ^q/xq 

r radial coordinate, m 

i’q circular duct radius, m 

t time, sec 

V axial steady flow velocity (function of radial coor- 
dinate), m/sec 

V Q axial steady flow velocity r free stream value, 
m/sec 

v acoustic particle velocity (vector), m/sec 

x axial coordinate, m 

y distance from the wall in the boundary layer, m 

y nondimensional distance from the wall in the 
boundary layer (y/6) 

Z acoustic impedance (P/vp, kg/m 2 / pec 

a complex radial eigenvalue (a = R 1 ^) 

~a dimensionless quantity (see Eq. (4-4), appendix) 

0 a(l-6/r 0 ) 


0 dimensionless quantity (see Eq. (A- 3)) 

6 boundary layer thickness, m 

d* boundary layer displacement thickness, m 

e dimensionless boundary layer thickness, 6/ r Q 

1 specific acoustic impedance (also used for optimum 

value with a finite boundary layer thickness) 

optimum specific acoustic impedance with zero 
boundary layer thickness (slip flow at wall) 

7) frequency parameter (fD/c) 

0 specific acoustic resistance, Real (£) 

optimum specific acoustic resistance when 6 = 0, 
Real (£ 0 ) 

0 boundary layer momentum thickness, m 

X sound wavelength, m 

v kinematic viscosity, m 2 /sec 

f nondimensional radial coordinate (r/r^) 

n 

p density, kg/m l 

<7 attenuation coefficient (see Eq. (10)) 

ctq maximum possible attenuation coefficient when 
6 = 0 

r propagation coefficient (see Eq. (10)) 

$ angular coordinate, rad 

(p phase of eigenvalue, a, deg 

X specific acoustic reactance, Imag (£) 

Xo optimum specific acoustic reactance when 6 = 0, 
Imag (£ 0 ) 

a? circular frequency, rad/ sec 

Subscripts; 

th 

j indicate j 1 radial mode 

r radial component 

Theoretical Model 

In the following sections the analytical approach will 
be outlined. The standard sound propagation theory will 
first be reviewed to establish the methods and terminol- 
ogy used in this paper. The method of obtaining the 
least attenuated mode optimum impedance will then be 
presented. Calculated results of the parametric study 
with boundary layers will next be presented followed by 
the correlation equations which unify all of the calculated 
results. Finally, some approximate expressions will be 
developed which are useful in illuminating the roles of 
the several variables considered here. 

The geometry and steady flow profile considered 
here are as shown in Fig. 1. The duct is circular with 
no splitter rings or hub. The boundary layer velocity 
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(5) 


profile is linear near the wall and has a 1/7^ power de- 
pendence outside of this linear region. The steady flow 
is assumed to be uniform in the central region of the 
duct outside of the boundary layer. 

Spinning Mode Propagation Theory With 
a Boundary Layer 

Certain assumptions are common in all of the theo- 
retical development which follows. Only in the boundary 
layer are steady flow gradients considered. The soft 
walls are assumed point reacting and of uniform imped- 
ance. These assumptions imply that no waves propagate 
within the liner backing cavity (prevented by honeycomb 
cells), and that any perforations are closely spaced com- 
pared to a wavelength of sound. No end reflections are 
considered; that is, the soft duct is sufficiently long to 
render the reflections at the duct termination unimpor- 
tant or the mode is well enough cut on so that no end re- 
flection will occur. A cut- on mode is defined here as a 
propagating mode; the frequency parameter ( 17 ) is suffi- 
ciently hi^i so that axial propagation occurs. Only very 
near cut-off will this assumption (no end reflection) be a 
problem. Below cut-off the modes carry very little en- 
ergy due to their out-of- phase pressure- v el ocity rela- 
tionship, and well above cut-off the end impedance^ is 
such that very little reflection will occur. , It Is further 
assumed that the sound pressure level is sufficiently low 
so that the linearized acoustic wave equation is valid. 

The linearized equations should be adequate except per- 
haps very near the rotor, when shock waves are gener- 
ated, or where near sonic flows exist. 

The wave equation with transverse (perpendicular 
to wall) steady flow gradients is given by: ^ 

JL i* + 2M + M 2 & - 2 p c™- ill - V 2 P = 0 
c 2 at 2 c ax at 0X 2 ar ax 


0 £ 


where in cylindrical coordinates 

v 2 p _0^P + 1_3P + J_ + 9^P _ 0 

dr 2 r 9 r j .2 9j£ 2 

Solutions are sought in the form: 

P = p(r)e i ^' im<i, - lkKx 
Combining Eqs, (1) to (3) yields 


(1) 


(2) 


(3) 


po;(l - MK) 8 r 

Defining a nondimensional radial variable 

i-s- 


where r Q - outer wall radius, Eq. (4) becomes 

£2+1 [1 , m dM l d p 

^2 $ L (1- MK) d|J d$ 


( 6 ) 


+ jim;) 2 [( 1 - MK) 2 - K 2 ] - ^ jp = 0 (7) 


where the frequency parameter 

n - - fD 

ire c 

Defining K as 

K - T - icr 

gives for the axial exponent, 

-ikKx = -k(cr + ir)x 


(S) 


(9) 


( 10 ) 


which indicates that, the attenuation is determined by cr, 
and the axial propagation speed is inversely proportional 

tO T. 

Uniform velocity region . In the central uniform 
velocity region ^ 1 - ( 6 /r Q ) or r ^ r Q - 5] , the Mach 
number gradient disappears (dM/d£ - 0) and Eq. (7) be- 
comes Bessels' differential equation with standard modal 
solutions 


p ] = V“f 


<u) 


where the j subscript indicates that an infinite set of 
modal solutions are possible. This discussion will be 
restricted to the least attenuated radial mode and the j 
will be dropped for brevity. The attenuation and propa- 
gation coefficients are related to the eigenvalue a by 


-iM f 


a + It — - 




1 -M 2 


( 12 ) 


££+±[n- 

J U ( 


2 Kr dM dp 
dr * i: l (1 - MK) dr J dr 

f 

, 2 


+ <K 


[U - MK ) 2 - K 2 




(4) 


where k = w/c, and use w r as made of the momentum 
equation 


The ratio of the pressure slope to the pressure 
evaluated at the boundary layer, is given by 


m - £ 


dp/U j 

p -U=l-(6/r 0 ) 1 -- 


J m<« 


(13) 
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which provides the boundary condition to start the Range- 
Kutta integration in the boundary layer and where 



profile provides a wall shear stress that is commensu- 
rate with the free stream velocity and boundary layer 
thickness assumptions. The velocity profile near the 
wall was derived from shear stress considerations given 
in Ref. 12. 


When the boundary layer is ignored (slip condition at 
the wall), continuity of particle displacement is used at 
the wall. The specific acoustic impedance result is 


where 


£ « 


G! 


W a > 

V“> 


m 


(15) 


Q = 1 + iM 0 (o- + ir) = 1 - MqK (16) 

Boundary layer region . When a velocity gradient 
exists, closed form solutions of Eq. (7) do not exist 
except for special oases such as a linear velocity pro- 
file. ^ Thus a Runge-Kutta integration of Eq. (7) was 
performed over the boundary layer with Eq. (13) as a 
starting point. More detail will be given in the Calcula- 
tion Procedure section. The boundary layer profile was 

— ■ 4 - 1 . 

assumed to be linear near the wall with a 1/7 1 power 
relation used over the rest of the boundary layer. 

At the wall M = 0, and the specific acoustic wall 
impedance can be calculated from 

t = zimni d7) 

dp/d^J^i 

where p and dp/di; at | - 1 are provided by the 
Kimgc- Kutta integration process. 

Boundary Layer Velocity Profile 


In the laminar sublaj r er near the wall the velocity 
profile can be expressed as 



where 

y * r 0 - r (19) 

is the distance from the wall, V Q is the free stream ve- 
locity, 6 is the entire boundary layer thickness, and v 
is the kinematic viscosity. 

Outside of the laminar region the velocity profile 
was assumed to be 



the velocity predicted by Eqs. (18) and (20) are equal and 
the transition from laminar to turbulent profile is made. 
A slight change in slope occurs at the transition point 
since no attempt was made to model the transition re- 
gion. 

The steady flow velocity and slope in the form nec- 
essary for use in Eq. (7) are in the laminar sublayer, 


It was concluded in Ref. 11 that a linear velocity 
profile with slip could be used to approximate the actual 
turbulent velocity profile. This would provide an ade- 
quate approximation if the boundary layer shape factor 
(6*/0) and displacement thickness (6*) were matched. 
This procedure greatly reduces the number of points 
required in the integration and is probably the sensible 
procedure to follow when a great number of calculations 
are to be made. In this paper, however, a limited num- 
ber of calculations were intended with the results being 
correlated for future use. The more conventional ap- 
proach using an adequate description of the entire bound- 
ary layer was thus used to eliminate the possibility of 
an unforeseen error due to an unnecessary assumption. 
This did, however, require that at least 100 points be 
used in the Runge- Kutta integration across the boundary 
layer. 

The following description of the boundary layer 


M = 


0,0225 M Q r 0 


(^f 


(l-«) 


(22) 


and 


dM . 


0.0225 M Q r 0 /M q c6\ 3//4 

~ \ v j 


(23) 


and in the turbulent region, 


M = 


= M 0 


( 1 "« ) 


1/7 


and 


(24) 


M 


dM 




(25) 


Note that in the laminar sublayer the boundary layer 
thickness (6) occurs in the Reynolds number (MqC6/ v ) 
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without duct size (Tq) scaling. Thus the laminar profile 
does not precisely scale with duct size and the size of 
the duct must be specified. For all calculations in this 
paper a duct radius of - 0 . 815 meter (36 in. ) was 
used. The scaling error should be small since the lami- 
nar region is only 2 percent or less of the boundary layer 
for typical Reynolds numbers of interest. 

It should be recognized that any velocity profile used 
is an approximation since the effect of an acoustic liner 
on the profile near the wall is not known at this time. 

Calculation Procedure 

The calculation procedure will be outlined first for 
the uniform flow case (no boundary layer) since the 
philosophy can be most easily explained and the boundary 
layer case is then a simple extension. 

Uniform Flow 

The intention here is to generate constant attenua- 
tion curves in the wall impedance plane for a particular 
mode. Where the closed curves of constant damping de- 
generate to a point, the maximum possible sound power 
attenuation is determined as well as the optimum wall 
impedance. This procedure for the uniform flow case 
without boundary layers is explained here and illustrated 
in Fig. 2 for M Q = -0. 4, 77 = 10, and m ~ 7. 

First a desired constant clamping value is selected 
(such as o - 0. 01 in Fig. 2). Next, a scries of propa- 
gation coefficients r values are selected. These are 
in turn (along with 0 ) used in Eq. ( 12 ) to calculate the 
eigenvalue a t then Q (Eq. (16)), and finally the acous- 
tic impedance £ from Eq. (15). Each impedance value 
(0 + ix - £) is normalized by 77 and plotted in Fig. 2, 

By selecting a new damping value (<r) the series of con- 
tours are generated. The optimum impedance and max- 
imum damping occur where the tear-drop shaped con- 
tours degenerate to a point. 

Boundary Layer 

A similar procedure to that of the uniform flow 
case is followed with a boundary layer, Again a and t 
are selected and Eq. (12) is used to calculate the eigen- 
value o'. Equations (13) and (14) are used to calculate 
l/p dp/d| at the boundary layer edge. This pressure- 
pressure slope value is used to start the Runge-Kutta 
integration through the boundary layer. This integra- 
tion transforms l/p dp/d£ at the boundary layer edge 
to l/p dp/d£ at the wall. Equation (17) is then used to 
calculate the wall impedance. This amounts to using 
the standard well known solution (Eq. (11)) in the uni- 
form flow region and deriving from this solution the 
boundary conditions necessary for the unknown Solution 


for propagation in the boundary’ layer described by Eq. 

(7). This procedure is very efficient since numerical 
integration is necessary over only a small portion of the 
duct (the boundary layer). Also modal correspondence 
or identification is quite simple since a predetermined 
mode is specified in the uniform flow region. This pro- 
cedure is similar to that used in Ref. 5. 

Equation (7) is solved by Runge-Kutta integration 
(using 100 steps) through the boundary layer using Eqs. 
(22) to (25) to Specify the Mach number and its slope. 
When the wall is reached in the integration, the pressure 
and its slope arc used in Eq. (17) to calculate the re- 
quired wall impedance. 

Sample results of this calculation are shown in 
Fig. 3. These are similar to the curves of Fig. 2 except 
that only a single damping value curve is shown for each 
boundary layer thickness. The optimum impedance is 
found within each one of these contours and can be esti- 
mated by the impedance value at the round tip of the tear- 
drop shaped contour. For the results presented in the 
next section the contours were further refined from those 
shown in Fig. 3. However, even for the sample curves 
shown here, it is evident that an extremely large drop in 
optimum resistance occurs when the boundary layer is 
increased. 

Results of Calculations 

The calculation procedure described in the previous 
section was exercised with a wide range of input varia- 
bles. Mach numbers from -0. 8 (inlet) to +0.6 (exhaust), 
spinning lobe numbers (m) from 3 to 50, frequency pa- 
rameters (tj) from 5 to 50, and boundary layers up to 
3(1 percent of the radius were used. The purpose of this 
section is just to present the results and point out the 
extreme magnitude of the boundary layer effect in the in- 
let (negative Mach number). The effect in the exhaust 
duct is much more moderate. Since the results are 
from a numerical integration where no closed form solu- 
tion exists, they may be looked upon as extremely re- 
peatable computer experiments of limited use, except to 
note trends, until they are unified and correlated. The 
correlation of the results will be performed in the next 
section. 

The results of interest in this paper are the optimum 
acoustic impedance and the maximum possible sound 
power attenuation obtained at this optimum impedance. 
Tills section contains results only for the least attenu- 
ated spinning modes. Higher order radial modes are 
deferred to a later section. Note also that the frequency 
parameter ( 77 ) is always selected so that the modes prop- 
agate. This was done to simplify the correlation pro- 
cedure and not because of limitations of the calculation 
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procedure. Also propagating modes are usually of most 
interest. 

The following results have been normalized by their 
counterpart with no boundary layer. This will increase 
the ease of use since the zero boundary layer thickness 
results are quite easily obtained. 

Inlet Conditions - Negative Mach Number 

Figure 4 shows the results for the optimum resis- 
tance ratio as a function of boundary layer- wavelength 
ratio for negative Mach numbers (inlets). This ratio 
compares the optimum resistance with and without a 
boundary layer, Note that the full boundary layer thick- 
ness is used here. If displacement thickness (6*) were 
desired instead it could be estimated as 6/8. Note the 
extreme reduction of optimum resistance with increasing 
boundary layer thickness particularly for large Mach 
numbers and frequency parameters. Note also that there 
is a frequency effect beyond that of wavelength which is 
accounted for in the abscissa (6/X). Also the Mach num- 
ber effect cannot be simply accounted for by using an 
effective wavelength [X(l + M^)l . 

Spinning lobe number (m) has an effect upon 0/ 
but this effect is not so dramatic as that of inlet Mach 
number or frequency. It will be found in the next section 
that these curves are all similar in shape with the roll- 
off value of 6/X being a function of the variables in- 
volved. 

A sample of the optimum reactance ratio is shown in 
Fig. 5 as a function of 6/X for m = 7, ri = 10, and for 
several inlet Mach numbers. The results (for a given 
Mq) increase at first to about double the 6 A = 0 value 
(X?Xq ~ 1) and then gradually fall off with increasing 
boundary layer thickness. Note that both x an d Xq are 
negative. 


boundary layer thickness, not much change occurs in the 
optimum resistance ratio (0.96 0/Bq — 1.05). Recall 
that for the same conditions (m - 7, u = 10) for an inlet 
(M 0 = -0.6) the optimum resistance ratio changes by 
several orders of magnitude as 6/X was varied over a 
similar range (see Fig. 4), Although the reactance ap- 
pears to change substantially, this is due to the method 
of presentation as a ratio. The actual reactance changes 
from about -0. 31 to zero in this boundary layer thickness 
range. The maximum damping ratio starts to increase 
somewhat at 6/X = 0. 25 (6/r Q = 0, 05) and has increased 
by only 30 percent at 5/x = 1 <6/r 0 = 0. 2) which repre- 
sents an extremely large boundary layer. 

It should be recognized that the model used here is 
not necessarily valid for an exhaust duct. Exhaust ducts 
are usually annular and the velocity profile is probably 
dominated by the rotor and stator properties rather than 
a simple boundary layer. The conclusion to be drawn 
here should be that all things being equal (except the flow 
direction) the exhaust duct optimum impedance is much 
less sensitive to sound refraction effects than is the inlet 
optimum impedance. 

Correlation of Results 

The correlating equations discussed below are de- 
rived in the appendix where assumptions leading to the 
correlations, deletions, and empirical corrections are 
discussed. The correlations are considered for inlets 
(negative Mach number) only since there the boundary 
layer effects are large while they are small for exhaust 
ducts. First, the complete correlation equations are 
presented which provide a very adequate unification of 
the calculations described in a previous section. Since 
these expressions are still quite complicated it is valu- 
able to look at more approximate expressions which bet- 
ter show the effects of the variables involved. 


The change in mode damping (o) is shown in Fig. 6. 
The maximum possible attenuation falls off very little 
until extremely large boundary layers are encountered. 
This implies that if the optimum impedance could be 
properly accounted for in a liner design, the peak atten- 
uation need not suffer much of a loss due to the boundary 
layer. 

Exhaust Conditions - Positive Mach Number 

Figures 7(a), (b), and (c) contain calculated results 
for optimum resistance, reactance, and damping for 
positive Mach number. These results were calculated 
with seven spinning lobes (m- 7), V - 10, and Mq- +0. 6 . 
The boundary layers contained in the span of the abscissa 
were from l/2 to 20 percent of the duct radius. 

In spite of die large Mach number and range of 


Complete Correlations 


The expression for the optimum acoustic impedance 
for spinning modes in a circular inlet (M^ ^ 0) with 
steady sheared flow is given by (see appendix) 


t - 


(1 + e)£ 0 


1 - 


in 


0 


(26) 


where is the optimum acoustic impedance without a 
boundary layer (slip flow), e ~ <5/rQ, and F is the ex- 
pression 


F _ tt(6/X) 

« 2< 1 + M 0> 2 


M 0 (1-2M 0 )-4(1 +2M 0 )(1 +M 0 ) 



(27) 


The quantity 


contains all of the boundary layer shape 
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information and for a l/7^ power velocity profile is ex- 
pressed as 


M„ 

/ M o\ 

1 + — 

i +-£ 

4 

\ 9 )■ 


(1 + jy 2 


(28) 


The quantity in the denominator of Eq. (27) is then 


2 

i 2 (i + m 0 ) 2 = i + -9 




(29) 


least attenuated spinning mode calculations, samples of 
which were previously shown in Fig. 6. This correlation 
is completely empirical in contrast to those for the re- 
sistance and reactance ratios. The damping correlation 
is shown in Fig. 10 by the solid curve while the previous- 
ly calculated points are shown as symbols. Again the 
results are shown as a ratio; maximum damping with a 
boundary layer divided by maximum damping with 5 = 0 
(slip flow). 

The correlation is given by 


The last term in the brackets in Eq. (27) is generally 
small (identically zero at = -0. 5) and most of the 
lobe number (m) effect is carried in the f ^ term in 
Eq. (26). 


a _ 1 

°0 5M?(6/\) 2 

1 + - 


(34) 


Equation (26) can be further reduced by taking real 
and imaginary parts yielding 


1 + e 

8 0 1 + 2Fx 0+ F 2 (o 2 + x 2 ) 


(30) 


The maximum attenuation ratio can also be calculated 
using Eq. (34), since the maximum attenuation is pro- 
portional to the damping coefficient for single modes. ^ 
Zero boundary layer thickness maximum attenuations for 
spinning modes are found in Ref. 1. 


JL 

*0 


(1 




X + 2Fx 0 + F 2 (o 2 + x^) 


(31) 


The zero boundary layer thickness input parameters 
to Eqs. (26), (30), and (31) can be estimated from Ref. 1 
as 



^ - 

(1 + M 0 ) 2 


(32) 


*0“ 


-: b -x!L- 

a + iy 2 


(33) 


valid for propagating modes. The quantities and 
are found in Table I (reproduced from Ref. 1 for 
convenience) for the least attenuated spinning mode with 
m lobes. 


Approximate Correlations 

It is desirable to simplify the resistance and reac- 
tance correlations and reduce them to more universal or 
general curves. Single correlating curves can be ob- 
tained for the optimum resistance and reactance ratios 
with only a small loss in accuracy. This will be done in 
two steps. First some approximations will be made to 
collapse the optimum resistance and reactance ratio 
correlations, and second some more gross approxima- 
tions will better reveal the effects of the major variables. 

Acoustic Liner Design Correlations 

In Table I it is apparent that « 4B^ which is 
fairly accurate except at m = 0 and 1 (the exact ratio of 
Bg/B^ could be used here). This implies that (see Eqs. 
(32) and (33)) 


d 0 ” - 4 *o 


(35) 


Equations (30) and (31) were used along with Eqs. 
(27) and (28) to calculate the normalized resistance and 
reactance ratios using all of the conditions shown on 
Figs. 4 and 5, The results are shown in Figs. 8 and 3 
as curves along with the previous more exact calcula- 
tions shown as point symbols. The agreement is seen to 
be very adequate for the resistance ratio and quite ac- 
ceptable for the reactance ratio. As discussed in the 
appendix, a slight empirical correction was made to the 
correlation equations to better fit the resistance ratio. 
No special measures were taken to fit the reactance. 


This can be used in Eqs. (30) and (31) to derive 


and 


e 0 1 - i FC 0 + (F0 O ) 2 


1 - 4F0, 


0 


x ° 1 + <*v 


(36) 


(37) 


Now it is apparent that the optimum resistance and re- 
actance ratios arc functions of only a single quantity, 
F0 q. This quantity is used as the independent variable 


A correlation was obtained for the damping of the 
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to plot Eqs. (36) and (37) in Figs. 11 and 12 which also 
contain the previously developed calculated points. The 
collapse of the resistance calculations onto a single 
curve is seen to be very good in Fig. ll, while some 
scatter is seen in the reactance ratio for very large 
boundary layers in Fig. 12. Again, this reactance error 
is not significant for liner design purposes. 

Note that F is negative for negative Mach number 
and thus -F was used in these two figures. The term 

(-4F0«) causes the initial increase in the reactance ratio 
u 2 

while the term (F0 q) causes the ultimate roll-off in 
both ratios. 

Equations (36) and (37) should be sufficiently accu- 
rate for inlet liner design purposes. The quantity F is 
calculated from Eq. (27) and is calculated from 
Eq, (32) using values from Table I. 


Further Approximations to Illustrate the Effect 
of the Variables 


It is interesting to look at the correlating equations 
for large (-F0 q). Equations (36) and (37) can be written 
in the limiting forms 


0 o <F0 O ) 2 

> 

*o FC, o J 


(38) 


(39) 


Thus 0/ varies with the inverse square power and 
x/Xq the inverse of (-F0 q). A comparison of Eqs, 
(38) and (39) with Figs. 11 and 12 show an adequate ap- 
proximation for large (-F0 q). 


A rough approximation for as seen in Table 1 
can be given by 


Q. 8 
1/2 


(m 4 0) 


(40) 


where m is the lobe number. This gives the approxi- 
mation (using Eq. (32)) 


0. 8 7] 


° m i/2 0 + M 0 ) 2 


(41) 


Further assume that the first term in Eq. (27) dominates 
the second term. This assumption is usually true. Also 
drop the empirically determined term l/2(l - 2 Mq), 
which is similar to one for most practical inlets. Then 
Eq. (27) yields 


F 


ttM 0 (6A) 
2^1 + M 0 ) 2 


(42) 


Next use Eq. (29) and ignore the M^/4 term. There- 
fore 

irM n 


0 / 5 


2 


Use Eq. (41) to obtain 


F0 0 * 


OAjtMqT, / fi 


m 1//2 (l + Mp) 2 ^ 


And finally from Eqs. (38) and (39) 


~ 0.6 m 


U +M 0 )^ 


and 


JL ~ - 


M 0 r,(5A) 
3.2-^m (1 +M 0 ) 2 


Xq MqJ] (6 A) 


(43) 


(44) 


(45) 


(46) 


Recall that these expressions are valid only for 
(-F0 q) > l, and for the least attenuated spinning mode. 
Some observations can be made concerning the effect of 
boundary layer refraction on the optimum impedance of 
acoustic liners. Note that reductions in these equations 
mean increased refraction effects. The extreme sensi- 
tivity to inlet Mach number and to the frequency parame- 
ter to a slightly smaller extent are evident in the equa- 
tions. This sensitivity was of course seen in the more 
exact calculated results. The higher the spinning lobe 
number the less the effect of sound refraction. 


It is also interesting that there is a frequency effect 
beyond that found in <5/\. Also it is not sufficient to use 
an effective wavelength = (1 + M^)X to account for 
the Mach number effect. 

To determine whether these approximate expres- 
sions (Eqs. (45) and (46)) can be used and in fact whether 
boundary layer refraction must be considered, the nomo- 
graph of Fig. 13 can be used. This nomograph deter- 
mines the values for which (-F0 q) - 1. For example, 
let Mq - -0. 5 and suppose a lobe pattern with 10 lobes 
is suspected. A value of - 2.25 is read off the 

ordinate. Now suppose that the boundary layer displace- 
ment thickness is 0,0028 r Q . The boundary layer thick- 
ness for purposes of the present correlation is 6/r rt - 

th v 

0.0224 (6 = 86* for 1/7 power velocity profile) . Thus 
tj 2 = 100 or 7] = 10. For a frequency parameter above 
10, the boundary layer refraction will cause sufficient 
change to alter the liner design and the regime of validity 
of Eqs, (45) and (46) is being approached. 


Refraction Effect on Higher Order Radial Modes 

The correlations presented in the previous section 
were for the least attenuated radial component of a spin- 
ning mode. It is of interest to determine whether the 
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correlation is sufficiently generalized to include propa- 
gating higher radial inodes. 

Equations (36) and (37) must be generalized to apply 
to the higher order radial modes. First Eq. (35) is al- 
tered to become 

^^0 ( 47 ) 

where ft was 4 for the least attenuated radial mode but 
is higher than 4 for the higher radial modes. The anal- 
ogies to Eqs. (36) and (37) can then be easily derived 
from Eqs, (30) and (31) and are expressed as 


fl o + 


and 


JL 

*0 


1 - RFe 0 

1 - — F0 O + <F0 O ) 2 
(ft U 


(49) 


These expressions are shown as curves on Figs. 14 
and 15 along with the exact calculations for the second 
and fifth radial modes of a seven-lobe spinning mode. 

For the second radial inode ft - 5. 9 and for the fifth 
radial mode ft = 10. 1 was used. The agreement be- 
tween the exact calculations and the correlations are 
seen to be excellent for the optimum resistance ratio 
over the entire (-F0 q) range, while the optimum reac- 
tance ratio deviates somewhat for higher values of 
(-F0 Q ), as did the optimum reactance ratio for the least 
attenuated mode. The reactance correlation could easily 
be improved if it were considered necessary. A slight 
error in reactance is not critical in a liner design since 
it will only slightly change the frequency at which peak 
attenuation occurs. Resistance is much more critical 
since it will determine if the maximum attenuation can 
be approached at all, and the effect of a boundary layer 
on the optimum resistance is much greater than on the 
optimum reactance. 


Because of the way in which the correlating equa- 
tions are normalized they seem to have general validity 
for any propagating mode. This normalization involves 
the use of the zero boundary layer or slip- flow solutions. 


Discussion of Boundary Layer Refraction Parameter 

It is apparent from the previous discussions of re- 
sults that the single parameter (-FOq) pretty well deter- 
mines the boundary layer refraction effects which be- 
come important when -F0 q > 1 for an inlet. The dis- 
cussion will be limited here to the alteration of the opti- 
mum resistance which is by far the dominant effect. 


-F0 n 


"-7tM 0 (5A) 

B e n 

_2I 2 (1 + M 0 ) 2 

(1 + M 0 ) 2 


(50) 


Increases in the bracketed quantities in Eq. (50) 
will increase the boundary layer refraction effect and the 
inverse is also true. The parameters which obviously 
increase refraction effects are increased boundary layer 
thickness (6), increased frequency (through bo th \ and 
tj), and increasingly negative Mach number (inlets). At 
high negative Mach number, the (1 + Mq) 2 quantity is by 
far the dominant term. Recall that roughly cancels 
one of the <1 + M Q ) 2 terms (see Eq. (29)). More subtle 
effects occur due to the term. The quantity is 
seen to decrease with increasing lobe number (m) in ta- 
ble I, and foi* any given value of (m), also de- 
creases with higher order radial mode number. Thus 
the effect of boundary layer refraction is diminished for 
higher order spinning and radial mode numbers. 

All of the boundary layer shape information is con- 
tained in the integral I ^ (see appendix). Since all of the 
results in this paper were for similar velocity profiles, 
no quantitative statements about lg will be made at this 
time. 

Although no claims have been made about the opti- 
mum resistance and reactance ratio correlations fitting 
the exhaust duct optimization calculations, the develop- 
ment in the appendix was general until the small empir- 
icism was added. Equation (50) does not contain the 
empiricism and could be easily altered for use with 
positive Mach number by dropping the minus signs. Re- 
call that only qualitative approximations are considered 
here. The quantities 6, X, r\ t and have the same 

effect as in the inlet case, and L still roughly cancels 

2 

one of the (1 + M Q ) terms. However, the remaining 
(1 + M q ) 2 term now decreases F0 O instead of increas- 
ing it as in the inlet case. For instance, the ratio of the 
F0 q values between minus and plus Mach numbers at 
IMq I = 0. 5 is a factor of 9. Thus the boundary layer in 
the exhaust duct would have to be nine times that in the 
inlet duct (all other terms being equal) to start produc- 
ing a significant refraction effect. The divergence be- 
tween the two cases (±Mq) is even more pronounced than 
the simple calculation above suggests. This can be seen 
in Eq. (36). For inlets, the -(l/2)F0 o term is actually 
positive and both terms in the denominator cause a fall- 
off in 0/0Q. However, for exhaust ducts -(1/2)F0 q is 
negative and first causes an increase in 0/ 0 q before 
the last term (FPq)^ takes over and causes a fall-off. 
This behavior can be seen in Fig. 7(a). 

Concluding Remarks 


It is interesting to look at the parameter in the 
following form 


With the optimum impedance information presented 
here and the spinning mode damping information of 
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Ref, 1, a design procedure for open circular inlet ducts 
was formulated which accounts for most of the important 
effects in the propagation of sound in a conventional 
acoustic liner. 

Some limitations of the procedure should, however, 
be recognized. Only propagating modes were considered 
in the development. When a mode is very near or below 
cut-off, some of the simplified expressions (such as 
Eqs. (32), and (33)) begin to lose validity. A rotation of 
the equal damping contours (see Ref. 1) occurs as well 
as the expansion <M Q < 0) or contraction (M Q > 0) con- 
sidered here. Also very near mode cut-off, the no end- 
reflection assumption may cause problems. A uniform 
boundary layer was assumed here. In practice, of 
course, the boundary layer is growing with length. 
Average values or stepped liners could be used to ap- 
proximately account for this growth. 

A 1/7^ power velocity profile away from the wall, 
merging into linear profile near the wall was assumed 
in the calculations. This was felt to be a reasonable ap- 
proximation to what may be encountered in a real inlet. 
All of the correlations are presented using the l/7 
power boundary layer thickness (5). If a boundary layer 
displacement thickness is available approximate 6 by 
5 = 88 *. 

Further study will be required to determine the 
sound refraction effects on the optimum impedance of 
exhaust ducts. Although very little sensitivity to refrac- 
tion was evident with the velocity profile used in this 
study, it is possible that a more realistic profile (for 
exhaust ducts) may produce a more significant effect. 

The calculations in this paper are for the optimum 
impedance values which are of use even in drawing in- 
ferences about attenuation results for off-optimum im- 
pedances. This can best be illustrated by an example. 
Suppose a calculation were being made using the param- 
eters given on Fig. 2 with a value of Q/t) = 1. 05 (which 
represents a resistance above the optimum) and x/tj ~ 

An attenuation coefficient <7 » 0.015 would be calculated 
if a very small boundary layer thickness were used. As 
boundary layer thickness is increased the optimum point 
moves down (to lower resistance, see Fig. 3) and the 
equal attenuation contours (cx constant) move down with 
the optimum. This causes lower damping contours to be 
swept past the selected impedance (1.05, 0) used in the 
calculation. Thus the calculated attenuation would mono- 
tonically decrease with increasing boundary layer thick- 
ness. Now suppose that a lower resistance had been 
selected for the attenuation calculation, for instance 
0/tj - 0. 55 (less than the optimum) and xA) = °* For 
very small boundary layers the calculation would yield 
a » 0. 01. Now as boundary layer thickness is increased 
the higher attenuation contours near the optimum would 


be swept by the impedance of interest (0 . 55, 0) and the 
attenuation calculation would at first increase. When 
the optimum impedance has passed below (due to in- 
creased boundary layer thickness) the selected imped- 
ance the attenuation calculation would then show a fall- 
off in <j. Thus the attenuation has first increased and 
then decreased as boundary layer thickness was in- 
creased. From these two examples it is seen that the 
trend of attenuation with boundary layer thickness for 
off- optimum impedances depends upon the relationship 
of tiie impedance to the optimum impedance. 

Conclusions 

A wide ranging parametric study was presented 
which considered the effect of boundary layer refraction 
on the optimum acoustic impedance of the least attenu- 
ated spinning mode and a limited number of higher radial 
modes in a circular duct. The following conclusions 
were made, 

1. Extremely large reductions in optimum acoustic 
resistance can occur due to the effect of an inlet steady 
flow boundary layer. The refraction effect of the bound- 
ary layer is less pronounced for the optimum acoustic 
reactance, 

2. Boundary layer refraction effects are increased 
in inlets for increasing Mach numbers, frequencies, 
and boundary layer thickness. Inlet refraction effects 
arc diminished for higher order circumferential and 
radial modes. 

3 . A correlation equation was derived, mainly from 
theoretical considerations, which adequately compares 
with the more exact calculations for an inlet, 

4. The correlation equation was found to be adequate 
for higher order radial modes, 

5. From the correlation, a boundary layer refrac- 
tion parameter was developed which indicates whether 
or not refraction should be considered in the design of a 
specific inlet suppressor. 

6. There is not much inherent loss of peak attenua- 
tion in an optimized design inlet due to boundary layer 
refraction. That is the maximum possible attenuation of 
a particular mode is a very weak function of boundary 
layer thickness until extremely high thicknesses are eiw 
countered. The optimum inlet design may have to be 
altered to account for refraction effects in order to 
maintain the peak sound attenuation at a high level. 

7. The effect of boundary layer refraction on the 
optimum, wall impedance in an exhaust duct is much less 
than for an inlet if all conditions (except flow direction) 
are the same. 
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Appendix - Development of Com- elation Equations 

The purpose of this development is to derive a cor- 
relation which can be used to unify the optimum imped- 
ance calculation results which were obtained by numeri- 
cal integration in the body of this paper. 

The starting point is an equation from Ref. 8 which 
expresses an approximate solution for the propagation of 
sound in a boundary layer velocity gradient and is valid 
for thin boundary layer thickness. This equation valid 
at the boundary layer edge is 


With 1/p dp/d£ a constant at the boundary layer 
edge, the right side of Eq. (A- 1) is used twice at zero 
thickness and at e to form the equality 

_ i7i7?A + € fel- + i?rrjA - (titj) 2 ] 

i7I7)A 0 _ =1 (A- 9) 

1 + it'^Al^ 

where is the optimum wall admittance at e = 0, 
while A is that at a finite e. The quantity (1- M Q K) 2 
was cancelled since K does not vary if e and t are 
nearly constant. Use was made of the expression 


i del 

-(1 - M 0 K) 2 [ikr Q A + £ ()JI 1 - «)] 

(A-l) 

P dlJj 

-=1- E l+lekr 0 AT 2 

where 


6 

€ = 

(A- 2) 


r o 

j8 = (kr 0 ) 2 K 2 + m 2 

(A-S) 


“ = _ikr o A + ( kr o> 2 

(A-4) 


1 = 1 

(A- 5) 


The quantity A is the acoustic admittance at the 
wall, and 1^ and are integrals over the boundary 
layer given by 



h = C (1 - M 0 KG) 2 <$ (A- 7) 

The coordinate y is the nondimensional distance from 
the wall (y/6). The function G represents the boundary 
layer shape and for a l/7 tl1 power low velocity profile is 
given by 

G(y) = y l/7 <0 sy =£ l) (A-8) 

Several simplifying assumptions and two fortuitous 
observations enabled the development of this complicated 
expression into a reasonably simple correlation expres- 
sion. First, it was noted that the damping (cr) and prop- 
agation (t) coefficients did not. significantly change at the 
least attenuation mode optimum impedance point as the 
boundary layer thickness was increased. This implies 
that the eigenvalue (note; o not T*) is also not a strong 
function of boundary layer thickness (see Eq. (12)) which 
in turn implies that the left hand side of Eq, (A- 1) can be 
considered constant for changes in e. This last impli- 
cation evolves from Eqs. (13) and (14). 


ki-Q = 70) (A- 10) 

Equation (A- 9) can be rearranged and using A = 1 /t 
there results 


(1 + e)£ 0 - ic? 


1 - 16717} 


h 


_ {TrnV_ 



(A- 11) 


Equation (A- 11) expresses the optimum wall impedance 
with a finite (but small) boundary layer thickness 
(e = in terms of the familiar (see Ref. 1) uniform 
or slip flow optimum wall impedance. This is useful 
provided the other terms in the expression can be eval- 
uated. 


The next step is to evaluate the integrals 1^ and I ^ 
which will require some simplifications. We will limit 
the development to propagating modes so that (see Ref. 1) 
a first order approximation for K is 

K = T -irjK i (A- 12) 

1 + M 0 

Thus I 1 and Ig can be written (using also (A- 8)) 

Ij = (j + Mq) 2 f & ( A- 13) 

Jo [l+M 0 (l-y l/7 )] 2 

and 

la = - 1 — - f 1 [l + ftyi - y l/7 )]“ dy (A- 14) 

(1 + M 0 > 2 - / 0 

To this author 1 s knowledge, cannot be evaluated in 
closed form but of course Ip can be. Upon integrating 
I-y numerically, the second fortunate occurrence was 
observed. It was found that up to |Mq I = 0.8 

Lb— (A- 15) 

3 % 

with an error of less than 9 percent. Equation (A- 15) 
was found to be approximately true also for a linear 
boundary layer with slip where both integrals can be 
evaluated in closed form. 
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The closed form expression for w *th a ve “ 
locity profile is 





(1 + M 0 ) 2 


(A- 16) 


Which is Eq. (28) in the main text. The integral ^ Is 
then taken as the reciprocal of Eq. (A- 16). 


Let F denote the quantity in the denominator of 
Eq. (A- 11) and use Eqs. (A-3) and (A-4) to yield 






(A- 17) 


2 

To evaluate F a second order approximation to K 
must be used (since m/irr\ also small). From Ref. 1 
one can derive (for propagating modes) 


K 


1 

1+m 0 


- (-^ L \ (cos 2(p + i sin 2<p) 
2 \mj/ 


(A- 18) 


where R and <p are the mode eigenvalue magnitude and 
phase at the optimum impedance. Since cp is usually 
small at the optimum (Table n of Ref. 1) 

cos 2(p > sin (A- 19) 


and the last term of Eq. (A- 18) will be ignored. Thus 


K 2 * 1 

(1 + M q ) 2 


f TJ \2 

— ) cos 2<p 


1 + M 0 


(A- 20) 


denominator (iek^Ay of Eq. <A-1) which may have im- 
plications in the use of Eversman 1 s result for moderate 
boundary layer thickness. 

The modified form of Eq. (A-ll) is then 


(1 + ^Q 

1 - 1FE 0 


(A- 23) 


which is Eq. (26) in the main text. Now if Eq. (A-23) is 
used with F from Eq. (A-21) retaining only the first 
term in brackets, a quite acceptable fit of the results in 
Fig. 4 is obtained. However, the agreement between the 
exact calculation and the correlation can be improved by 
using small empirical corrections. This is not sur- 
prising due to all of die approximations made to this 
point. The simplest possible corrections were made 
(linear corrections) and somewhat more sophisticated 
curve- fitting procedures could probably improve the 
correlation even more. 

A Mach number multiplier was inserted into the 
first term of Eq, (A-21), This correction term was 
l/2(l - 2 Mq) which ranges from 0. 7 to 1. 3 as Mach 
number varies between -0. 2 and -0. 8. 

The second empirical correction was made to the 
second term of Eq. (A-21) which involves the quantity 
cos 2 (p which is the real part of the square of the 
eigenvalue. This is a function of the circumferential 
lobe number as well as the radial mode number. Rather 
than inserting the actual values of R and cp t the term 
was considered empirically as a function of the lobe 
number (m) . The empirical substitute obtained was 
(l/2)m 2 . 


and Eq. (A- 17) can be written 

M rt \ / _ \ 2 


P = 2tt(oA) 

h 


1+ T/ ( ; 005 2tp / \2 
' 9 / + W /mA 

,) w 


4(1 + M 0 r 


(1 - M 0 ; 


With these two empiricisms inserted, Eq. (A-21) 
becomes (also use mg/9 « 1) 


(A-21) 


F - ■ 


jeML 


4 ¥ 1 + M o> 


M 0 (l - 2M„) - 4(1 +2M 0 )(1 +M 0 >(^y 


where 


was used. 

This is about as far 
be carried (Eq. (A-21». 
must be made with the more exact calculations made in 
the main text. When the real part of Eq. (A- 11) (using 
F in the denominator) was compared to the calculated 
points in Fig, 4, it was immediately evident that the 
second term in the numerator of Eq. (A- 11) had to be 
deleted. This term (-ie^y causes 0/0^ to go nega- 
tive at moderate boundary layer thickness. Deleting 
this term in effect is deleting the small term in the 


- 2 f $ 
tj U 


(A- 22) 


as the approximate theory could 
From here on compari sons 


(A- 24) 

which is Eq, (27) in the main text. 

The intent in this appendix was to develop a bound- 
ary layer refraction effect correlation which was as 
theoretically based as possible, hi this way the major 
effects of the parameters would appear in their correct 
functional form with a minimum of guesswork which 
usually accompanies a correlation. Small empirical 
corrections were then applied to improve the correla- 
tion. The final form of the equation would probably 
never have been derived from purely empirical curve- 
fitting. 
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TABLE 1. - OPTIMUM RESISTANCE AND REACTANCE 
COEFFICIENTS FOR SPINNING MODES 
IN CIRCULAR DUCTS 


m 

0 

1 

3 

r> 

7 

10 

20 

50 


a. 890 

0. 644 

0.448 

0. 357 

0. 303 

0. 251 

1 

0. 170 

0.097 

R 

X| 

0.38 

0.221 

0. 121 

0. 095 

0. 076 

O.OOO 1 

0.038 

0.020 
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Figure 1. - Geometry and steady-flow velocity profile. 
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Figure 2. - Constant attenuation and propagation contours 
for the least attenuated spinning mode using displacement 
boundary condition. Inlet Mach number, M 0 = -0.4, lobe 
number, m=7; frequency parameter, n « 10; boundary 
layer thickness, 6/r 0 * 0. 
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Figure 3. - Constant attenuation contours for several 
boundary layer thicknesses (6/r 0 ) for the least attenu- 
ated spinning mode. 
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Figure 4. - Effect of boundary layer on inlet optimum resistance 
ratio for the least attenuated spinning mode with lobe number, 
frequency parameter, and Mach number as parameters. 
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BOUNDARY-LAYER THICKNESS TO WAVELENGTH RATIO, 6/A 


Figure 5. - Effect of boundary layer on inlet optimum reactance 
ratio for the least attenuated spinning mode with Mach number 
as a parameter. 



BOUNDARY LAYER THICKNESS TO WAVELENGTH RATIO, 6/A 

Figure 6. - Effect of boundary layer on inlet optimum damping 
ratio for the least attenuated spinning mode. 
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(c) DAMPING AT THE OPTIMUM. 


Figure 7. - Exhaust duct optimum resistance, reactance, and 
damping as a function of boundary- layer thickness. Least 
attenuated mode, lobe number m = 7, frequency parameter 
tj = 10, Mach number M 0 = 0.6. 
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Figure 8. - Optimum resistance ratio, comparison between exact 
calculations and correlation for inlets for the least attenuated 
spinning mode. 
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BOUNDARY- LAYER THICKNESS TO WAVELENGTH RATIO, 6/X 

Figure 9. - Optimum reactance ratio, comparison between exact 
calculations and correlation for inlets for the least attenuated 
spinning mode. 



Figure 10. - Optimum damping ratio, comparison between 
exact calculation and correlation for inlets for the least 
attenuated spinning mode. 




Figure 11. - Inlet optimum resistance ratio for the least attenuated 
spinning mode, unification of calculations using boundary- layer 
refraction parameter. 
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Figure 12. - Inlet optimum reactance ratio for the least attenuated 
spinning mode, unification of calculations using boundary- layer 
refraction parameter. 



Figure 13. - Nomograph for onset of significant boundary- layer 
sound refraction effect on optimum resistance. 
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Figure 14. - The effect of inlet boundary layer on the optimum resist- 

z ance ratio for higher order radial modes. 
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Figure 15. - The effect of inlet boundary layer on the optimum re- 
actance ratio for higher order radial modes. 



